On the Aizenman exponent in critical percolation 
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The probabilities of clusters spanning a hypercube of dimensions two to seven along one axis of 
a percolation system under criticality were investigated numerically. We used a modified Hoshen- 
Kopelman algorithm combined with Grassberger's "go with the winner" strategy for the site per- 
colation. We carried out a finite-size analysis of the data and found that the probabilities confirm 
the Aizenman's proposal of the multiplicity exponent for dimensions three to five. A crossover to 
the mean-field behavior around the upper critical dimension is also discussed. 



Percolation occurs in many natural processes, from 
electrical conduction in disordered matter to oil extrac- 
tion from field. In the latter, the coefficient of oil extrac- 
tion from oil sands (the ratio of the actually extracted 
to the estimated oil) can be as much as 0.7 for light oil 
and as low as 0.05 for viscous heavy oil. An increase 
in this coefficient by any new point requires appreciable 
investment. Additional knowledge about the percolation 
model could reduce the amount of additional investment. 

A remarkable breakthrough in the theory of critical 
percolation was established in the last decade thanks to 
a combination of mathematical proofs, exact solutions, 
and large-scale numerical simulations. Recently, Aizen- 
man has proposed a new exponent that describs the prob- 
ability P(k, r) of a critical percolation (i-dimensional sys- 
tem with the aspect ratio r being spanned by at least k 
clusters |BJ, 
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hxP(k, r) oc —ay fc^ r, 



(1) 



where ad is a universal coefficient depending only on the 
universality class, and £ = d/(d — 1). 

In two dimensions, Aizenman's proposal (|l|) was 
proved mathematically jfl, confirmed numerically [Q, 
and derived exactly || using conformal field theory and 
Coulomb gas arguments. This exponent seems to be re- 
lated to the exponents of two-dimensional copolymers j| . 
In three dimensions, proposal (Q) was checked numeri- 
cally in H and, more recently and more precisely, in 

The upper critical dimension of percolation is d c = 
6, which follows from the comparison of the exponents 
derived on the Cayley tree with those satisfying scaling 
laws (see, e.g., JtJ and ||). The fractal dimension Df of 
percolating critical clusters is equal to 4 above d c , and 
the number of percolating clusters becomes infinite for 
d > d c . This fact would imply that ( = at d = 6 if 
we supposed (rather naively) that Aizenman's formula 
applies at the upper critical dimension. Supposing that 
this is true and taking into account that the values of £ 
for d = 2 and d — 3 are, respectively, 2 and 1.5, we can 
place all three points on the straight line ( = (6 — d)/2, as 
depicted in Fig. [l| We can then estimate the respective 
values of ( for d = 4 and d = 5 to be ( = 1 and £ = 0.5; 



FIG. 1: Variation of Aizenman exponent £ with the space 
dimension d as predicted by Aizenman (circles and dotted 
line), claimed by Sen (dashed line) and discussed in the text 
(solid line). 



these values are far from those predicted by Aizenman's 
formula giving 4/3 and 5/4, respectively. In contrast, 
based on simulations, Sen || claims that £ = 2 for all 
dimensions from two to five. 

The main purpose of our simulations is to estimate 
the exponents for the dimensions from two to six with an 
accuracy sufficient for distinguishing between the values 
predicted for d — 4 and d — 5 by the Aizenman's formula 
and a naive application of cluster fractal-dimension ar- 
guments and by the straight-line fit, as discussed above. 

In the rest of the paper, we briefly summarize the high- 
lights of our study, then present some details of our re- 
search, and finally discuss the results for the Aizenman 
exponent and the physics of a crossover from the Aizen- 
man picture to the mean-field picture. 

Our main results can be summarized as follows: 

1. Modified combination of the Hoshen-Kopelman al- 
gorithm and Grassberger's strategy. We use the Hoshen- 
Kopelman (HK) algorithm JTcJ to generate clusters and 
Grassberger's "go with the winner" strategy to track 
spanning clusters. We add a new tag array in the HK 
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algorithm, which allows the reduction of the tag memory 
order from L d to L^ 1 , where L is the linear size of the 
hypercubic lattice. As a result, the amount of memory is 
about two orders less for large values of L, and the pro- 
gram is about four times faster — the complexity of the 
algorithm is compensated by the lower memory capacity 
needed for for swapping to and from the auxiliary array. 

2. Efficient realization of combined shift-register ran- 
dom number generators. We use an exclusive-or (©) com- 
bination z n of two shift registers: 

%n = ^n-9689 © ^n-5502) 

Vn = 2M-4423 © 2M-2325, z n = x n © y n (2) 

(see pi and the references therein) . We reduce the com- 
putational time for generating random numbers by a fac- 
tor 3.5 through an efficient technical modification: we use 
the SSE command set that is available on processors of 
the Intel and AMD series starting from the Intel Pentium 
III and AMD Athlon XP. 

3. Extraction of the exponents for dimensions three 
to five. We first use finite-size analysis to estimate the 
logarithm of the probability P(k, r) in the limit of infi- 
nite lattice size L. We then fit data as a function of the 
number of spanning clusters k to obtain the Aizenman 
exponent 

4. Confirmation of Aizenman's proposal. The esti- 
mates of the exponent C for the dimensions d = 2,3,4, 
and 5 coincide well with those proposed by Aizenman. 

5. Qualitative interpretation of Aizenman's conjecture. 
Cardy interpreted Aizenman's result qualitatively in two 
dimensions on the basis of assumption that the main 
mechanism for reducing the number of percolation clus- 
ters is that some of them terminate. The same result can 
be derived for the cluster confluence (or merging) mech- 
anism. This means that, in low dimensions, the percola- 
tion clusters consist of a number of closed paths (loops), 
while, in higher dimensions, clusters are more similar to 
trees. Indeed, it is well known that the probability of 
obtaining loop becomes lower for higher dimensions and 
goes to zero in the limit of infinite dimensions (Cayley 
tree) § §. 

6. Crossover to mean-field behavior. We found evi- 
dence that the probability of clusters spanning a hyper- 
cubic lattice tends to unity in the limit of high dimen- 
sions, as it follows from the well-accepted picture. We did 
not find any dramatic changes in the probabilities around 
the upper critical dimension d c = 6, but rather found evi- 
dence for a crossover. Therefore, Aizenman's formula (|l|) 
can also apply to dimensions higher (but not too much 
higher) than the upper critical dimension and describe 
approximately the probabilities of spanning clusters in 
large, though finite-size systems. 

We follow with the details of the critical percolation, 
simulations, and data analysis. 
Spanning probability. 
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FIG. 2: The logarithms of the probabilities of exactly k clus- 
ters P(k,r;L) (+) and of at least k clusters P+(k,r; L) (x) 
for the dimension d — 4 and the number of clusters k — 5 as 
functions of the aspect ratio r. The linear size of the hyper- 
rectangle is L = 16. The solid line is the linear approximation 
to \nP(k,r; L) on the interval r — [1.5; 5.0]. 

We can define the probability P(k, r; L) that k clusters 
traverse a <i-dimensional hyperrectangle [0, x [0, Lr] 

in the Lr direction 0] . Provided that the scaling limit ex- 
ists (this was proved recently by Smirnov for the percola- 
tion in plane the probability P(k, r) can be defined 
as the limit of P(k, r; L) as L — > oo. Aizenman proposed 
that P(k,r) should behave according to (|l]) in dimen- 
sions from three to five. The validity of formula (11]) for 
the percolation in plane was well established in pL |3|. 

Numerical results (|| and ||) for the exponent ( for 
critical percolation on cubic lattices seems to confirm 
Aizenman's proposal for the value of £ — 1.5 

Actually, we could consider the probability P(k,r) as 
the probability of obtaining k clusters at the distance r 
from the left side of the hyperrectangle if clusters grow 
to the right. Only two processes can change the number 
of clusters: cluster merging and cluster terminating. 

The differential dP of the probability is 

dP oc P(k,r) kW-^k dr, (3) 

where the right-hand side represents the product of the 
probability P(k, r) and the differential of the total border 
hyperarea of k clusters, each with the hyper area differ- 
ential fcVtd-i) dr. This expression follows from the fact 
that area unit of measure is proportional to the charac- 
teristic transverse length of "infinite" clusters. There- 
fore, the transverse area remains constant as k changes, 
while the longitudinal length increment in these units is 
oc k^-l^^dr. Integrating (||), we recover probability (Q). 
Thus, P(k, r) describes the probability that k clusters do 
not merge together. 

The same probability could be obtained by the process 
of cluster termination, as given by Cardy in plane [0, 
which can easily be extended to dimensions d > 2. 
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d k Lrnin Lmax SL Vc 

2 1-5 16 256 16-32 0.59274621(13) 

3 1-6 8 64 4,8 0.3116080(4) 

4 1-6 8 48-56 4,8 0.196889(3) 

5 1-6 4 32,24 4,8 0.14081(1) 

6 1-6 4 15-16 3-5 0.109017(2) 

7 1-4 4 10 1 0.0889511(9) 



Ref. 

IP I 
[L4| 

[15| 

[L5| 

[16| 

[16| 



TABLE I: Minimal L m in and maximal L max linear sizes of the 
percolation lattice and the interval SL between two consecu- 
tive values of L depending on the dimension d and number of 
clusters k. The values of p c are taken from the references in 
the last column. 
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This means that the exponent £ cannot be larger than 
the one proposed by Aizenman, and £ = d/(d — 1) is the 
upper bound for the exponent. 

Algorithms and realizations. 
The classical realization of the HK algorithm jl(| requires 
memory for two major structures: an array for keeping 
a (d— l)-dimensional cluster slice and a tag array. The 
total memory required by the algorithm is oc L d ~ x + 
p c rL d , where p c is the site percolation threshold value. 
Therefore, for large rL, one of the main advantages of the 
HK algorithm (i.e., relatively low memory consumption) 
is negated by the second term. Our modification of the 
original algorithm allows the memory for the tag array 
to be reduced to about 3p c L x . 

Instead of keeping all tags in memory and selecting new 
tags with increasing tag numbers, we create two arrays, of 
which one keeps the tag value and the other one keeps the 
number N of the slice where the corresponding tag was 
last used. When we build a cluster, we update this array 
with N = iV curr e n t for the tags used. If N < N culrcnt — 1, 
then this tag is not on the front surface of the sample, 
and it will never be used again so that we can, therefore, 
reuse it. We note that cluster size information should be 
taken into account before reusing the associated tag, if 
the size information is required. 

We use the "go with the winner" strategy as fol- 
lows. If the system has k spanning clusters for some 
aspect ratio r — nSr, it is stored in memory and is 
grown for Sr. If the resulting configuration has k span- 
ning clusters, it is stored, and the growth process con- 
tinues. Otherwise, we return to the previously saved 
state. Using this procedure, we calculate the probabil- 
ity Pi(Sr) that the system propagates at the distance 
Sr from the position r = (i — 1)6 r. Finally, we ob- 
tain P(r — nSr) — Y\a=i Pi{$ r )- By choosing sufficiently 
small values of Sr, we can achieve rather high probabili- 
ties of Pi{5r) (which can be determined from a few real- 
izations), while the total probability may be very small 
(down to oc 10~ 100 in our case). 

The random number generator was optimized for the 
SSE instruction set as follows. Because the length of all 



FIG. 3: Plot of s(k;l) for k = 3 clusters in the dimension 
four as a function of 1/(1/ + Lo) with the fitting parameter 
Lo = 4.12 (□) and for k = 4 clusters in the dimension six 
with Lo = 4.03 (o). Straight lines result from the fitting to 
the corresponding data as discussed in the text. 



k 






d 






3 


4 


5 


6 


7 


1 


-1.377(1) 


-1.774(3) 


-1.859(9) 


-1.76(2) 


-1.48(4) 


2 


-6.919(6) 


-6.330(15) 


-5.57(6) 


-4.73(8) 


-3.55(11) 


3 


-13.655(15) 


-11.64(4) 


-9.95(12) 


-8.27(12) 


-6.25(16) 


4 


-21.47(3) 


-17.77(6) 


-14.65(20) 


-11.95(25) 


-9.3(3) 


5 


-30.23(3) 


-24.02(8) 


-19.9(3) 


-15.75(30) 




6 


-40.02(6) 


-31.0(1) 


-25.0(3) 


-22.7(2) 





TABLE II: Values of s(k) for different numbers of clusters k 
and dimensions d for site percolation on hypercubic lattices 
with periodic boundary conditions in directions perpendicular 
to the spanning direction. 



four RNG legs is {a^j^ij,} > 4, the nth step of the RNG 
does not intersect with the (n+3)th step. Therefore, we 
can pack four consecutive 32-bit values of {x^^i^}} 
and {y n -{ a \b }} into 128-bit XMM registers, process 
them simultaneously (see Eq. (|2|)), and thus obtain z n , 
2 n +i, z n+ 2, and z n +3 within one RNG cycle. 
Data analysis. 

The lattice size was varied from L m i n to L max with the 
step SL. In Table [j], particular values of the simulation 
parameters are presented together with the interval of 
the number of clusters k depending on the dimension d. 
The direct result of the simulations is the probabilities 
P(/c,r;i) that exactly k clusters connect two opposite 
surfaces (separated by the distance rL) of the rectangle 
with size L d_1 in the "perpendicular" direction in which 
we apply periodic boundary conditions. We use the val- 
ues of the site percolation thresholds on hypercubic lat- 
tices from |l^]-[^6), as shown in Table |. 

Data analysis consists of three steps. First, we com- 
pute the slope s(L) of lnP(fc, r; L) for a given dimension 
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k this paper exact from [3] from [5 



1 -0.6541(5) -0.6544985 -0.65448(5) 

2 -7.855(3) -7.85390 -7.852(1) 

3 -18.32(1) -18.3260 -18.11(15) 

4 -32.99(3) -32.9867 

5 -51.83(2) -51.8363 

TABLE III: Values of s(k) in two dimensions for different k 
calculated in thisnaper, using the exact Cardy formula 
and estimated in ]5[ for site percolation on a tube. 



d, number of clusters k, and linear lattice size L. An ex- 
ample of such a function is given in Fig. for In P(5, r; 16) 
in the dimension four. We also plot the logarithm of the 
probability P + (k,r;L) — J2k'>k r \ L) of the event 
that, at least, k clusters span the (hyper)rectangle at the 
distance rL. To calculate s(L), we use data only in the 
interval of the aspect ratio r between 1.5 and 5. We note 
that the probability of five clusters spanning a rectangle 
with the linear size L = 16 at the distance 5 • 16 = 80 is 
extremely small sa 10 -52 . 

Second, we compute probabilities in the limit of an infi- 
nite system size L, fitting slopes s(k) with the expression 
(see Fig. §) 



s(k;L) = s(k) + 



D 



(L + L y 



(4) 



where B, t, and Lq are fitting parameters || |l7| H - 
The resulting values of the slopes s(k) are presented in 
Table ||. The number of runs used to compute each 
particular entry in Table [□] varied from 10 6 to several 
tens for higher dimensions. 

We checked the accuracy of our simulations, as well as 
the validity of the approach in general for site percola- 
tion on a square lattice. Table |lj shows a comparison 
of our results for the slope s with the exact values and 
with earlier simulations, in which the other modification 
of the HK algorithm, but not the Grassberger strategy, 
was used. We note that our results coincide well with the 
exact results and give a higher accuracy for larger values 
of k in comparison with the previous numerical results, 
despite the smaller computation time used. Our data for 
k = 1 is less accurate because of the smaller statistics 
(10 6 runs, compared to 10 8 samples in ||). This is a di- 
rect demonstration of the effectiveness of the Grassberger 
strategy for large values of k. 

Finally, we use values in Table || to determine the 
Aizenman exponent £ by fitting data in each column to 



s = A(k 2 

in two and three dimensions, 
berger |6), and to 



k ) p/2 



(5) 



as proposed by Grass- 



d A 


ko 


P 


2 2.090(4) 


0.244(5) 


2.0012(10) 


2.0940(5) 0.2489(7) 2 


3 2.81(4) 


0.64(4) 


1.489(7) 


2.757(2) 


0.587(3) 


3/2 


4 3.06(20) 


0.41(6) 


1.315(30) 


2.949(5) 


0.373(3) 


4/3 


5 2.8(1) 


0.40(4) 


1.24(3) 


2.78(2) 


0.38(2) 


5/4 


6 2.8(8) 


0.5(3) 


1.12(14) 


2.41(5) 


0.33(6) 


6/5 


7 1.4(10) 


0.08(116) 1.4(4) 


2.03(12) 


0.50(13) 


7/6 



TABLE IV: Values of the fitting parameters A and ko and 
the power p as defined in Eqs. (g) and (^) for the dimension 
d. 



in higher dimensions. Here, A, ko, and p are fitting pa- 
rameters. We take only the leading behavior in k into 
account. 

Spanning, proliferation, and crossover to mean-field be- 
havior. 

The results of the final fit to (§) and (|) are shown in 
Table 



A(kP - k ) 



(6) 



: IV . The second row for each particular dimension d 
is the fit with the power p fixed to the Aizenman exponent 
value. This is done to check the fit stability. Indeed, the 
values of A and kg coincide within one standard deviation 
for the dimensions two to five. 

The larger deviations of parameters for the dimensions 
six and seven may be attributed to appearance of cluster 
proliferation — the number of clusters is known Q to grow 
as L d ~ 6 in dimensions d > d c — 6. We plot the coefficient 
aa (defined by expression (Q)) in Fig. | as a function of 
the dimension d. The probability of exactly one cluster 
spanning at the given distance r becomes smaller as the 
dimension increases from two to five and larger for larger 
dimensions, as can be seen from the first row (k = 1) of 
Table |l| and from the lower curve dependence in Fig. ^|. 
For any fixed d, the value of ad approaches some limit 
for the dimensions two to five and k > 2, which suggests 
the value of the corrections to the leading behavior in k 
(see Eqs. (g) and (§)). 

The fact that the value of (, which we formally ex- 
tracted from our data for d = 6, more or less coincides 
with £ = dj (d— 1) = 6/5, as formally calculated using the 
Aizenman expression, may be interpreted as an indica- 
tion that the number of clusters depends logarithmically 
on the lattice size L. One can expect that the logarith- 
mic behavior is visible only for somewhat larger values of 
L than we have used so far (see Table |). With the val- 
ues of L of the order we have used in simulations, we see 
effectively the same picture as for the lower dimensions — 
clusters spann according to the Aizenman formula. This 




FIG. 4: The coefficient ad (as a function of the dimension d) 
extracted from the probabilities P(k, r) for different numbers 
of clusters k. 



means that at small (or moderate) values of L, the main 
mechanism is as discussed above: cluster merging and 
terminating. And only at sufficiently large system sizes 
we will see cluster proliferation. An indication of that 
can be seen from the values of ad in the dimension seven 
in Fig. The probabilities become closer, and this can 
be attributed to cluster proliferation and treated as a 
crossover to the mean-field behavior. 
Discussion. 

The results have shown the validity of the Aizenman's 
proposal in the dimensions three to five (results on 
plane were already proved rigorously) and do not sup- 
port Parongama Sen claims based on their simulations 
(Fig. 0) . We have found evidence for cluster proliferation 
for the dimension seven. The analysis can be extended to 
the number of spanning clusters to distinguish exponen- 
tial decay with the system size of the number of clusters 
for the dimension five, logarithmic growth of them for the 
dimension six, and linear growth for the dimension seven. 
The same technique can be used to establish numerically 
such a crossover to the mean-field picture, although a 
significantly longer computational time, than we used, is 
needed for this. In fact, the linear growth of the multi- 
plicity of spanning clusters for seven-dimensional critical 
percolation was confirmed numerically in preprint pl| 
posted at arXiv preprint library a few days after our 
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